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UhiTRAl 1 


Tlir mt lucnco i’l combined st ites i>t st ri"- 1 ' on the sheat response 
il*'U)' nat.rlil prtinlp.il direi t inn-, lit ott-ixis. unidiroi t uni 1 1 , i nn’ii’s 1 1 e 
. oupons is exuninod lor ststons whose nitrix nalerfil obexs the von Mises 
' i old condition. Sin ti in.il v s I s is ni’tiv.itid bx trends in exper ineut .i 1 
dil i observed tor it Ieist l vo composite sxstcras th.it mdii ite dex iat ton 
trori pure stie.ir behavior ilong mterl.il d i rei t ions tor <• n lous ott-txi.s 
i out t eiir.it ions. 

Tin* vield timet ion tiir plane stress ot .1 t r.msverse I x isotropii 
i oripos tie 1 imiita . onsist me ot stilt, linoarlv elastic t ibers uid 1 
YOU Mises matrix miteri.il is tornxilated in terns of Hill's elastii stress 
i oik out i at ion tutors and i single plastii lonstriint par met or . The above 
ire subsequent 1 \ evaluated on the basis ot observed ive rage lamina and 
constituent ri spouse tor the Aveo SS05 boron-epoxx system. It is shown 
that Im lusion ot residuil stresses in the vield tunetion together with the 
itn orporat ion of IHibex and lllllier's concept ot general wed vield stress 
foi anisotropii nedi i nt the constitutive equation correctlv predii ts the* 
trends observed in experiments. The incorporat ion ot the strong axial 
stress interaction t'eci'ss irx to predict the correct trends in the shear 
response is directlx trued to the high residual axial stresses in the 
mtrix induced during fabric ition of the composite. 
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INTRODUCTION 


The use of the off-axis tensile specimen for determination and/or 
verification of elastic properties and fracture strength of unidirectional 
composite materials is well documented, eq. Sih [1], As most composite 
materials exhibit varying degrees of nonlinear response prior to failure 
in an off-axis tension test, this specimen would also appear well suited 
for verification of the various assumptions in a number of nonlinear 
theories recently proposed In the literature, as well as some of their 
consequences. An indirect outgrowth of certain assumptions employed in 
one such theory, for instance, is a recently proposed method for determi- 
nation of nonlinear lamina shear stress-strain response on the basis of 
the 10* off-axis tension test (2], 

The nonlinear response of unidirectional composites is a complex 
phenomenon which can be caused by a number of diverse mechanisms such as 
nonlinear constituent behavior, damage accumulation through fiber or 
matrix cracking or debonding, fiber rotation or any combination of the 
above to mention the better known causes. The extension of results of 
the uniaxial test along principal material directions to constitutive 
models for irultlaxial loading situations is not always straightforward 
due to the possibility of various stress interactions taking place at the 

microlevel. Foy [3] for example demonstrated the influence of the trans- 
it 

verse normal stress component on the shear strain response (and vice 
versa) along the macerial axes for composites with stiff, linearly elastic 
fibers and ductile matrices using a finite element, micro-mechanical anal- 
ysis. He showed that for combined normal shear loading in a fixed ratio 
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yielding initiated sooner and the extent of inelastic behavior was more 
pronounced than for either normal or shear stress acting alone (Fig. 1). 

' A nunber of nacroaechanical constitutive nodels have been proposed in 
the literature 14-81 to describe the nonlinear response of lanina with 
stiff, elastic fibers and significantly nore coapliant and often ductile 
matrices. The early attenpts [4,5] disregarded the possibility of stress 
interaction in the nonlinear range during combined loadings, while in the 
later formulations this effect was taken into account, often in an 
heuristic manner. One such model due to Hashln et al (6] is reniniscent 
of the total deformation theory of plasticity in that the stress interac- 
tion is accounted for through a loading function, based upon an invariant 
formulation, which is used to determine the different functional forms of 
the various strain components. The loading function is limited to quad- 
ratic terms in transverse normal and shear stresses hus eliminating the 
possibility of the axial normal stress influence on the transverse and 

shear strain behavior. 

Data generated by Cole and Pipes 19] are employed in this paper to 
re-examine and study the extent of stress interaction in off-axis unidi- 
rectional boron-epoxy coupons with the aid of Hashin's model and a micro- 
mechanical approach utilizing a yield function concept based on Hill's 
method [10] for determining the elastic stress concentration factors. 

The results of the micromechanical analysis are subsequently employed in 
Dubey and Hillier’s generalized yield stress formulation [11] to study 
the nonlinear shear stress-strain response along the principal material 
directions in off-axis tension coupons. Several fiber orientations are 
considered including the 10* coupon which has been proposed for the 
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STRAIN 


STRESS-STRAIN CURVES FOR UNIMRECTI0NA1. BORON/EPOXY (50% KIBFR 
VOLUME FRACTION) UNDFR COMBINED LOADS. NORMAL :SHFAR LOAD RATIO 
IS 8:3, REF. 3. 
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Jctcrmin.it ion of pure shear response. 

■Hie Influence of stress interaction on the shear response in the 
off-axis tensile tests of Cole and Pipes is shown in Fig. 2 along with the 
pure shear response obtained from unidirectional tubes. These experimental 
results show that the extent of inelastic behavior is not a monotonic 
function of fiber orient' . ion. The inelastic shear strains (for a speci- 
fied shear stress) decrease in magnitude as the fiber angle is increased 
from 15* to 30°, but as the fiber angle is increased further to 45* the 
inelastic strains increase. (Similar behavior has been observed by the 
authors for the graphite-polyimide off-axis coupons.) One of the goals 
of this paper is to provide an explanation for such behavior. 

In the first part of the paper the model of Hashin et al is briefly 
outlined. It is then used to predict the shear stress-strain response 
along material axes using unidirectional data taken from Ref. (9) . An 
empirical observation is made that modification of Hashin's model by 
incorporation of an axial normal stress interaction tern significantly 
improves correlation with the experimental data. Based on the above 
results, micromechanical formulation of the yield function for the 
composite in the presence of stiff fibers is employed along with Dubey's 
generalized yield stress concept to examine the effect of the stress 
component on yielding in the matrix phase and its consequences on the 
nonlinear shear response along material directions. It is shown that 
incorporation of the presence of residual curing stresses in the composite 
leads to the correct prediction of trends in the experimental shear stress- 
strain data of Cole and Pipes. 
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■ PURE SHEAR DATA (REF 9) 
A 30° COUPON DATA (REF 9) 
a45° COUPON DATA (REF. 9) 
• 15° COUPON DATA (SANDHU) 
0 15° COUPON DATA (REF. 9) 


0 -005 .010 .015 .020 .025 

ENGR. SHEAR STRAIN 

FIG. 2. EXPERIMENTAL SHEAR STRESS-STRAIN RESPONSE ALONG MATERIAL PRINCIPAL 

«n»y T «c?rM 0R VARI0US OFF-AXIS TENSION COUPONS. AVCO 5505 BORON- 
truAY oYoTEM* 




STRESS INTERACTION MODEL 


toll 0W||||> i J.isslcnl concepts of deformation theory of plasticity 
Hashin et al [6] assume that the inelastic part of strain can be expressed 
In the following runner* 

p _ > 

C ij " °ij £ ij (h) (no summation) (1) 

where L is some general quadratic loading function of the stresses. 

Since composites which employ stiff fibers such as boron or graphite 
behave linearly (or nearly so) under loads applied in the fiber direction, 
the plastic response is assumed to be limited to transverse and shear 
strains. Furthermore, the effect of the stress in the fiber direction 
on the remaining strain components is neglected and thus Eq. (1) is 
expressible as follows: 



g 12 , ,2 , n 2 . 
2C 12 f 12 (a °22 + 0O 12 ) 



E^ f 22 (a S 2 + eg2 12> 


( 2 ) 


where the interaction term f ’ 12 a 22 bas been excluded on the basis of 
material symmetry arguments and only strains relevant to laminate analysis 
retained. The functional form of ^2^) an< * ^12^ is determine ^ from 
uniaxial tension and torsion tests along material directions. It is noted 
that no expli.J. mention is made of the yield function and Its relationship 
to the loading function in this formulatior despite apparent separation of 
the total strain into elastic and inelastic or plastic portions. Combining 
plastic and elastic portions of strain, the relevant total strain components 
become : 
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q n v i2 

11 E n' E n 22 


t^o 11 + 
t ll 11 


c - -l 2 + 
12 2G 12 


..•i-B(ft)’-ft)'] 

4 [ft 


where the constants M, N, c^, and are obtained fron the Ranberg-Osgood 
approximation of uniaxial tests. 
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TH EO R ET I C AL- EX P E R I MENTAL-C ORRELATTON 


Data presented by Cole and Pipes on the Avco 5505 boron-epoxy system > 
liave been employed to determine the Ramberg-Osgood parameters in Hashin's 
model. These are given in the appendix. Predicted shear response in the 
material coordinate system for the various off-axis fiber orientations is 
plotted in Fig. 3. While Hashin's model predicts nonotonically decreasing 
response with increasing off-axis angle, the experimental data points to a 
reversal around 30* off-axis angle (Fig. 2). 

The concept of stress interaction in an elastoplastlc material can be 
illustrated by considering the shear stress-strain response along axes 
inclined at some angle with respect to the loading direction of an isotropic 
tensile specimen. The hypothetical isotropic material is taken to have the 
same shear properties as the considered boron-epoxy system for illustration 
purposes. Using the isotropic theory (in this case the flow and defor- 
mation approaches are coincident) the shear response along the inclined 
axes is given by; 

°12 . °12 r 2 . 1 2 n . _ 2„,, _ 2~ / °12,N-1 

C 12 " 2C^ + 2cq [ 3 + 3 (tan 0 + cotan 0)1 ( “ } (4) 

with yield governed by: 

°12 I j + J (tan 2 0 + cotan 2 0)) ” S 2 , O<0<9O® (5) 

/3 

(or alternatively, o 12 “ — j S sin20) where S is the yield stress in 
in pure shear. Eq. (4), plotted in Fig. 4, predicts monotonically 


Increasing shear response along rotated axes with increasing angle between 
0* and 45’ with reversal occuring at 45®. This follows from the symmetry 
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of the yield function (Eq. 5) about the 45° angle. For comparison purposes 
llashln’s theorv (Fig. }) yields the following relation for the shear res- 
ponse along material axes inclined at the same angle as the reference axes 
In the hypothetical isotropic .specimen 


°12 ° 


+ _ni [1 + ( V tan 2 0 j 2 ( liV-l 


N-l 


^ *'^’12 ^^12 


(b) 


y y 

The absence of the resolved shear stress reversal in the above model can be 
directly traced to the absence of the axial stress (o^) Interaction term. 
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MICROKECHANICAL CONSIDERA 1 -OHS 


Assuming that the matrix can be modelled as a von Hlses solid that 
controls the nonlinear response of the composite, the yield function of 
the lamina can be determined from the knowledge of the stress concentration 
factors at yield. We take the view that these stress concentration factors 
are related to the average elastic concentration factors through a single 
plastic constraint parameter. This appears to be partially substantiated 
by various finite-element, microraechanical studies, c.f. Dvorak and Rao 
(12], which have revealed that although yielding usually starts at the 
fiber/matrix interface, the stress vector remains nearly radial under 
proportional loading. The yield zone on the other hand spreads rapidly 
throughout the entire matrix with Increasing external tractions. 

The elastic stress concentration factors can be determined from the 
knowledge of the average lamina as well as constituent response following . 
Hill's outline [10]. The relationship between average matrix and external 
stresses at yield, including residual stresses, Is then assumed to be 
expressible as follows: 


B 11 B 12 


B 21 B 22 


where subscript m refers to matrix, r residual and no subscript external 


stresses. Also, where 


-, and the elastic 


stress concentration factors B.., ,B,. are given in the appendix. 


Here, K p is the previously mentioned plastic constraint factor to be 
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determined from experimental data, and the single scale factor A is a 
consequence of transverse isotropy of the composite. 

Employing the above formulation, von Mises yield function for the 
transversely isotropic lamina in plane stress becomes: 


3 (B ll~ B ll B 21 +B 21 )o ll K p + 3 (B 12 _B 12 B 22 +B 22 )o 22 K p 


+ 3 ( 2 B 11 B 12 + 2B 21 B 22 -B 11 B 22 -B 21 B 12 )a u o 22 K p 

+ y (i - A )o llmr (2B ll" B 21 )o ll K p + 3 (1 “ A )o llmr (2B 12~ B 22 )o 22 K f 

. _2 -2 2 v 2 1 ... 1 .- .2 

+ B 66°12 K p ” K “ 3 [(1 A^°llnr 


( 8 ) 


where K is the yield shear stress of the matrix material. 

The importance of the residual stresses cannot be overemphasized. 
Employing a single inclusion composite cylinder model and temperature 
dependent material properties for this particular system, the average cur- 
ing stresses have been estimated to be: 



Although time dependent response has not been considered due to insuf- 
ficient data, the reduction in residual stresses due to viscoelastic 
response will be partially off-set by the use of a single inclusion com- 
posite cylinder modelwhlch generally underestimates induced stresses. 

The above yield function (Eq. 8) is plotted in Fig. 5 in terras of the 
resolved shear stress along the material axes and the off-axis angle. 


NORMALIZED SHEAR YIELD STRESS 
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FIG. 5. YIELO SHEAR STRESS OF AVCO BORON-EPOXY SYSTEM ALONG' THE 

MATERIAL PRINCIPAL DIRECTIONS OF AH OFF-AXIS TENSION COUPON' 
AS A FUNCTION OF THE FIBER ORIENTATION. 
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Also plotted are the yield functions of an equivalent isotropic material 
and a transversely isotropic lamina with no residual stresses. It is 
clearly seen that the reversal of the shear stress-ptrain response evident 
in Cole's and Pipes' data can be explained on the basis of initial yielding 
in the matrix. To do this we formulate a generalized yield stress function 
according to Dubey's and Hillier's suggestion as follows: 

5 “ c ij°ij + ^iVk7ij°w 

and assume that the plastic shear strain is a function of 0 . The coeffi- 
cients are determined from the yield function given by Eq. (8) 

and the pure shear response. These are given in the appendix. The result- 
ing shear stress-strain curves along material axes for the various off-axis 
angles are presented in Fig. 6 indicating the extent of stress interaction. 
It is seen that the use of the 10° off-axis specimen for determination of 
pure shear response can lead to noticeable deviations for a material system 
that can be modelled by the assumed set of governing equations. 
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FIG. 6. PREDICTED SHEAR STRESS-STRAIN RESPONSE ALONG MATERIAL PRINCIPAL 

DIRECTIONS FOR VARIOUS OFF-AXIS TENSION COUPONS BASED ON DUBEY I 

AND HILLIER'S GENERALIZED YIELD STRESS CONCEPT AND MICRO- \ 

MECHANICAL CONSIDERATIONS INCLUDING RESIDUAL STRESSES. AVCO 
5505 BORON-EPOXY. SYSTEM. ! 
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CONCLUS IONS 


The devlationb from pure shear response along material axes in an off- 
axes tension test for the Avco 5505 boron-epoxy system have been explained 
to a large extent on the basis of nicromechanical formulation of a vield 
surface end a generalized yield stress concept. The model is applicablc- 
to systems with stiff, lineirl\ elastic fibers and ductile matrices whose 
response can be approximated by that of a von Mises solid. The trends 
evident in the data of Cole and Pipes for the shear response along the 
material principal directions are traced to the strong axial stress (o^) 
interaction which in turn follows from relatively significant residual 
stresses. The model predicts that in the absence of these stresses the 
Influence of axial stress on nonlinear shear stress-strain response for 
low off-axis fiber orientations is only significant in the 0®-5° off-axis 
range. When the residual stresses are included, however, this influence 
is noticeable up to approximately 15®. This is interesting in view of the 
recently proposed test method for determination of pure nonlinear shear 
stress-strain response on basis of the 10° off-axis tension test. 
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APPENDIX 


i) Rambert, ’ t tnneterb for the Avco 5505 Boron-Epoxy system 


M 

0 (|>bi) 

y 

N 

t (psi) 

1.434 

26,155.2 

4.430 

9,593.9 


li) Elastic stress concentration factors B,,,....B,, 

11 oo 


B„ =- 
11 v 

m 


(E.-E ) 2 -(v E -v E ) 2 
f m r m n i 


12 v 


(v f E n~ V m E f )(V f l: ir V 12 E f ) 

,: n 


<E f- E J 2 - l ”fV”„ E f 2 > 


( ”fV”o E f HE (- E 22> 


21 v 


22 


' E 


22 v 
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66- ' v 


(E f- E n> - (V f E m-'’n E f ) ' 


(E f" E n )(v f E ir V l2 E f ) 
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(E -E ) 2 -(v E -v E f ) 2 
i m i m n t 


(E f - E m )(E f - E 22 ) 


it 


rx2n)( c ±) 

V=r c i ' 


( E f-E n )(Ef-Eu) 

E 11 


(E f~ E m )(v f E ll~ V 12 E f ) 

E 11 


- (v f E m - v m E f )(E r E n> 

E 11 


~ (v f E n- V n E f )(v f E ir V l2 E f ) 

• E ll 
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* 0.19 x 10 psi 

= 0.49 > 10 6 psi 

v = 0.31 
m 


= 24.0 x 10 psi 
Ej. = 58.0 x 10^ psi 


v = 0.50 

in 


0 12 = 0.88 X 10 psi 
E^^ = 30.1 x 10^ psi 
E » 2.87 x 10 6 psi 


v 12 - 0.225 


in) Determination of generalized yield stress. From experiment 

e p - °12 r q ! 2. N ~ 1 - 

12 2G 12 ( T } * A1S ° ° ° /2D 66 °12 ln con tracted notation. Thus 

. Jli ( _ V 5 ) N - X 

12 !g 12 t /!dT7 

y 66 

From longitudinal tension and compression we obtain: = £ ,/d — i n 

contracted notation, where 

e . 5 lV 

<vv 

X t’ x c being the nagnitudes of yield stress in tension and compression 
respectively. Similarly in transverse tension' and compression 

C Z “ n • 


( V Y t> 

(Y c + V 
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From the lnv.iri.nu e of «' ue also have: 


/JO 

JfT~ - bb 

ni ( X t ) (1+0 


/2D 

-tfTT- <?-> — - 

t (l+r,) 

where S is the yield stress In pure shear. Thus the function — 


becomes; 


0 1 ,S . f - , ,S . 1) - 

, SET ’v S ' u V < ,+, » ” 

y 66 


V -2 -2 

S * °11 S ‘ °' f 2 -2 

(fo — ^ + —^7 + ^ - 

t (l+O* T (l+n) 


The parameters S, X [t , t" and n have been evaluated from the viold 
function Riven bv F.q. (8) and tlie experimental data. The plastic 


constraint factor has been determined on the basis of welding (i.e 
deviation from linearit>) in transverse tension. 





